Combining first KamLAND results with solar neutrino data 
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We consider the impact of the recent KamLAND data on neutrino oscillations, the first terrestrial 
neutrino experiment that can probe the solar neutrino anomaly. By combining the first 145.1 days 
of KamLAND data with the full sample of latest solar neutrino data we find an enhanced rejection 
against non-LMA oscillations, allowed only at more than 4a with respect to LMA. Furthermore, the 
new data have a strong impact in narrowing down the allowed range of Am^ inside the LMA region. 
In contrast, our global analysis indicates that the new data have little impact on the location of the 
best fit point. In particular the solar neutrino mixing remains significantly non-maximal {3a). 



I. INTRODUCTION 

In a recent paper the first results of the KamLAND 
collaboration became public Q . These data contain pre- 
cious information on the neutrino oscillation hypothe- 
sis which has been advocated to account for a num- 
ber of neutrino experiments involving solar and at- 
mospheric neutrinos and which indicate that neutri- 
nos are massive and that neutrino flavor mixing oc- 
curs HHIilllllllllSllSlllllllll. The Kam- 
LAND experiment is a reactor neutrino experiment with 
its detector located at the Kamiokande site. Most of 
the Pe flux incident at KamLAND comes from plants 
at distances of 80 — 350 km from the detector, mak- 
ing the average baseline of about 180 kilometers, long 
enough to provide a sensitive probe of the LMA solution 
of the solar neutrino problem [13. ITsl The Kam- 

LAND collaboration has for the first time measured the 
disappearance of neutrinos traveling to a detector from 
a power reactor. They observe a strong evidence for the 
disappearance of neutrinos during their flight over such 
distances, giving the first terrestrial confirmation of the 
solar neutrino anomaly and also establishing the oscil- 
lation hypothesis with man-produced neutrinos. More- 
over the parameters that describe this disappearance in 
terms of the oscillations of the electron neutrino type to 
another, are consistent with latest pre-KamLAND deter- 
minations [Ulii El iS m m lill of solar neutrino 
parameters. 

In this note we analyze the implications of these fun- 
damental results by combining the KamLAND data with 
data from solar neutrino experiments. We will assume 
CPT conservation and for simplicity we consider a two- 
flavor massive neutrino oscillation framework. In Sec. |n] 
we analyze the impact of the KamLAND results by in- 
cluding the full information on the spectral distribution 
of the observed events. Subsequently, in Sec. Illll we per- 
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form a global fit that combines the full KamLAND and 
Chooz reactor data sample |24| with the full solar neu- 
trino data as included in Ref. [13 ■ In Sec. IIVI we check 
the stability of the results with respect to changes in the 
statistical analysis, and we summarize in Sec. 



II. SIMULATION AND ANALYSIS OF 
KAMLAND DATA 

In KamLAND the target for the flux consists of a 
spherical transparent balloon filled with 1000 tons of non- 
doped liquid scintillator. The anti-neutrinos are detected 
via the inverse neutron /3-decay 

De + P ^ + n . (1) 

In Fig. 5 of Ref. the spectral data are given in 13 
bins of prompt energy above 2.6 MeV. To simulate the 
KamLAND data we calculate the expected number of 
events in each bin for given oscillation parameters as 

r (2) 

J2 <l>j{E,)Pj{E,,Am^, 0) / dEe R{Ee,E'^) . 

Here R{Ee,E'^) is the energy resolution function and 
Ee,E'^ are the observed and the true positron en- 
ergy, re spectively , and we use an energy resolution of 
7.5%/v/£^(MeV) _1]. The neutrino energy is related to 
the positron energy by E,^ — E'^ + A, where A is the 
neutron-proton mass difference. The integration interval 
over Eg is determined by the prompt energy interval in 
each bin. The neutrino spectrum (l){E^) from nuclear re- 
actors is well known, we are usin g th e phenomenological 
parameterization given in Refs. 25, 26]. We adopt the 
average fuel composition for the nuclear reactors given in 
Ref. 01 . Note that possible effects due to time variations 
in the fuel composition have been shown to be small ^2^ . 
The sum over j in Eq. (0) runs over 16 nuclear plants, 
taking into account the different distances from the de- 
tector and the power output of each reactor (see Table 3 
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Figure 1; Allowed regions at 90%, 95%, 99% and 99.73% 
C.L. (2 d.o.f.) and best fit point from KamLAND spectral 
data. The colored regions and the star are obtained by using 
the of Eq. Q based on the Gaussian approximation, the 
hollow regions and the dot correspond to the of Eq. Ulljl 
implied by Poisson-distributed data. 



Figure 2: Allowed regions at 90%, 95%, 99% and 99.73% C.L. 
(2 d.o.f.) from the combined analysis of solar, Chooz and 
KamLAND data. The hollow lines are the allowed regions 
from solar and Chooz data alone. The star (dot) is the best 
fit point from the combined (solar+Chooz only) analysis. 



of Ref. [23) • The relevant detection cross section <t{E^) 
is given in Ref. ^2§\. In the 2-ncutrino framework the 
survival probabihty for the neutrinos coming from the 
reactor j is given by 

Pj{E,,Am^e) = l-sm^2esm^^^^. (3) 

4-C/j/ 

The normahzation factor / in Eq. |(5Jl is determined in 
such a way that for the case of no oscillations we obtain 
a total number of events of 86.8, as expected from the 
Monte-Carlo simulation used in Ref. 

For the statistical analysis we use the x^-function 

= ^(iVf - N°'^^)S-/{N^'^ - iVf ^) . (4) 

The observed number of events N°^^ in each bin can be 
read off from Fig. 5 of Ref. 0. In the covariance ma- 
trix S we include the experimental error in each bin ai 
(obtained from the same figure), which we assume to be 
uncorrclated, and the systematic error dsyst = 0.0642 jlj 
implied by the uncertainty on the total number of events 
expected for no oscillations: 

S,,^afS,,+al,,N^'^N;^. (5) 

This x^ definition assumes Gaussian distribution of the 
data. For the discussion of an alternative analysis based 
on poisson distributed data see Sec. IIVI 



III. RESULTS AND DISCUSSION 

Our results are summarized in Figs. ^ [21 and |S1 In 
Fig. n we show the allowed regions of the oscillation pa- 
rameters obtained from our re-analysis of the KamLAND 
data. It is in good agreement with the analysis performed 
by the KamLAND group, shown in Fig. 6 of Ref. 
This gives us confidence on our simulation of the Kam- 
LAND data and therefore encourages us to use it in a 
full analysis combining also with the solar data sample. 
Figs. 121 and O show the corresponding results obtained in 
a combined fit of the full KamLAND data sample with 
the global sample of solar neutrino data, as well as the 
Chooz result. The solar data we are using and the details 
of our solar neutrino analysis are given in Ref. |l7j . 

First of all, we have quantified the rejection of non- 
LMA solutions and found that it is now more robust. For 
example, for the LOW solution we have Axlow-lma ~ 
26.9, which for 2 d.o.f. (Ato^ and 9) corresponds to a 
relative probability of 1.4 x 10~^, assuming Gaussian er- 
rors. A similar result is also found for the VAC solution. 
Apart from selecting out LMA as the unique solution of 
the solar neutrino problem we find, however, that the 
new reactor results have little impact on the location of 
the best fit point: 

tan^ e = 0.46, Am^ 6.9 x 10~^ eV^ (6) 

In particular the solar neutrino mixing remains signifi- 
cantly non-maximal, a point which is rather important 
for model-building. Indeed bi-maximal mixing models 
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Figure 3; Ax^ versus Am^ and tan^ 6. The red dashed hue refers to KamLAND alone. The green dot-dashed hne corresponds 
to the full reactor data sample, including both KamLAND and Chooz. The blue solid line refers to the global analysis of the 
complete solar and reactor data. 

are disfavored poj while models where the solar mixing with a Ax^ = 4.5 with respect to the best fit point given 
can be non-maximal [30l | are preferred, as before. This is in Eq. (|SJ|. This ambiguity might be resolved when more 
not in contradiction with the fact that KamLAND data KamLAND data have been collected (see e.g. Refs. 
alone prefer maximal mixing P| , since such preference Isil l3^ ) . 
has no statistical significance. Indeed, one can see from 
the right panel in Fig. O that Ax^ is rather flat with re- 
spect to the mixing angle for tan^ > 0.4. This explains 
why the addition of the KamLAND data has no impact 
whatsoever in the determination of the solar neutrino os- 
cillation mixing. The allowed 3a region we find for 9 is: 

0.29 < tan^ 9 < 0.86, (7) 

practically identical to the pre-KamLAND range given 
in Eq. (4) of Ref. ll7j]. 

On the other hand, the new data do have a strong im- 
pact in narrowing down the allowed range of Am^. From 
the left panel of Fig. O one can read off that KamLAND 
data alone provides the bound Am^ > 8 x 10^^ eV^, 
whereas the CHOOZ experiment gives Am^ < 10"'^ eV^, 
both at 3tT. Hence, global reactor neutrino data pro- 
vide a robust allowed interval for Am^, based only on 
terrestrial experiments using artificial neutrino sources. 
However, combining this information from reactors with 
the solar neutrino data leads to a significant reduction of 
the allowed range: As clearly visible in Fig. |21 the origi- 
nal LMA region is now split into two sub-regions. From 
Fig. Owe obtain at 3a (1 d.o.f.) 

5.1 X 10"-^ cV^ < Am^ < 9.7 x lO""^ eV^ (8) 

1.2 X 10"* eV^ < Am^ < 1.9 x 10"* eV^ (9) 
The local minimum in the region ^ occurs for 

tan^ 9 = 0.42, Am^ = 1.4 x 10~* eV^ (10) 



IV. STABILITY OF THE STATISTICAL 
ANALYSIS 

The current KamLAND data sample consists of 54 
anti-neutrino events, which are distributed over the 13 
energy bins. This leads to rather small numbers of events 
in each bin. The 5 bins with highest energies contain no 
event at all. In such a case the use of a x^-function based 
on Poisson statistics might be appropriate. In order to 
check the stability of our results we have performed also 
an analysis by using [3^ 



aN}^ - Nf'' + N°^' In ' " 



1-a 



•'syst 
(11) 

where the term containing the logarithm is absent in bins 
with no events. We minimize with respect to a in order 
to take into account the overall uncertainty of the theo- 
retical predictions. 

The analysis of KamLAND data using Eq. (|11() is 
shown in Fig. ^ as the hollow lines. We observe that 
this analysis is somewhat less constraining compared to 
the analysis based on the Gaussian x^ of Eq. @. One no- 
tices that smaller values of the mixing angle are allowed, 
especially at high convidence level. Let us note, however, 
that the allowed regions from the Gaussian analysis are 
in better agreement with the analysis done by the Kam- 
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Figure 4: As Fig. |21 but using for the analysis of KamLAND 
data the of Eq. Ull^ implied by Poisson distributed data 

LAND group. This is the reason why we prefer to use 
this method for analyzing KamLAND data. The better 
agreement with the original KamLAND analysis might 
be related to the fact that the inclusion of the information 
on the experimental errors provided by the KamLAND 
collaboration in Fig. 5 of Ref. l| can only be included by 
means of a Gaussian x^-function, as in Eq. In this 
way it is possible to take into account the asymmetric 
errors and the error bars in bins where the number of 
events is zero. 

However, we note that the determination of Am^ is 
rather stable, only the constraint on the mixing angle is 
somehow affected. Since the bound on the mixing an- 
gle in the combined analysis is dominated by solar data, 
we expect the difference between the two methods to be 
small after combining KamLAND with solar data. The 
results of this exercise are shown in Fig. 01 Comparing 
this figure with Fig. [21 we find indeed, that the result 
is very similar. The location of the best fit point and 
the 90%, 95% and 99% C.L. regions around the best fit 
point are nearly identical. However, the local minimum 
does not appear at the 90% C.L., though its location is, 
again, very stable. Some small differences are visible for 
the 99.73% C.L. contour. 

To summarize, although there are some notable dif- 
ferences between the allowed regions obtained by assum- 
ing Gaussian or Poisson x^-functions for the KamLAND 



data taken alone, the differences are very small when 
combined with solar data. This illustrates the robust- 
ness of our results against variations in the statistical 
analysis. 

V. CONCLUSIONS 

We find that among all previous oscillation solutions 
to the solar neutrino anomaly, the new reactor results 
from the KamLAND experiment single out the LMA so- 
lution, rejecting all other oscillation solutions at a signif- 
icant level. Furthermore, we find that already these first 
145.1 days of KamLAND data lead to a significant im- 
provement in narrowing down the allowed range of Am^ 
when combined with solar neutrino data. The original 
LMA region now is split into two relatively narrow is- 
lands around the values of Am^ = 6.9 x 10~^ eV^ (best 
fit point) and Am^ = 1.4 x 10~* eV^ (local minimum). 
However, our full analysis indicates that the new data 
have little impact on the determination of the mixing 
angle. In particular the solar neutrino mixing remains 
significantly non-maximal (Str). 

Before closing, let us note that we have considered 
here only the simplest case of two neutrinos. Analyz- 
ing in detail the impact of the KamLAND results on 
three-neutrino oscillation scenarios ,3^] and the resulting 
constraints is beyond the scope of this short note. The 
improved determination of Am^ can also play an interest- 
ing role in probing fine details of solar physics i mat- 
ter effects l^^; probing electro-magnetic neutrino prop- 
erties ^3 (see also "s^) or testing CPT invariance in the 
neutrino sector ^39] . Similarly, the nailing down of LMA 
as the solution has also implications for non-oscillation 
solutions to the neutrino anomaly |40j | in terms of spin- 
flavor precession l4l|. |43. non-standard interactions 
or neutrino decay |44| . Clearly none can now be the 
leading explanation to the solar neutrino anomaly [45l| . 
although a detailed evaluation must be performed to de- 
cide, in each case, to what extent these solutions are now 
rejected. 
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